Unified ab initio approach to bound and unbound states: 
no-core shell model with continuum and its application to ^He 
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We introduce a unified approach to nuclear bound and continuum states based on the coupling of 
the no-core shell model (NCSM), a bound-state technique, with the no-core shell model/resonating 
group method (NCSM/RGM), a nuclear scattering technique. This new ab initio method, no-core 
shell model with continuum (NCSMC), leads to convergence properties superior to either NCSM or 
NCSM/RGM while providing a balanced approach to different classes of states. In the NCSMC, the 
ansatz for the many-nucleon wave function includes: i) a square-integrable A-nucleon component 
expanded in a complete harmonic oscillator basis; ii) a binary-cluster component with asymptotic 
boundary conditions that can properly describe weakly-bound states, resonances and scattering; 
and, in principle, Hi) a three-cluster component suitable for the description of, e.g., Borromean halo 
nuclei and reactions with final three-body states. The Schrodinger equation is transformed into a 
system of coupled-channel integral-differential equations that we solve using a modified microscopic 
R-matrix formalism within a Lagrange mesh basis. We demonstrate the usefulness of the approach 
by investigating the unbound ^He nucleus. 

PACS numbers: 21.60.Dc,24.10.Cn,25.10.-fs,27.20.-fn 
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I. INTRODUCTION 



One of the central goals of nuclear physics is to come 
to a basic understanding of the structure and dynamics 
of nuclei, quantum many-body systems exhibiting bound 
states, unbound resonances, and scattering states, all of 
which can be strongly coupled. Ab initio (i.e., from 
first principles) approaches attempt to achieve such a 
goal for light nuclei. Over the past fifteen years, ef- 
ficient techniques such as the Green's function Monte 
Carlo (GFMC) ah initio NCSM j^]. Coupled Cluster 
Method (CCM) [1^ or nuclear lattice effective field the- 
ory (EFT) |6l| have greatly advanced our understanding 
of bound-state properties of light nuclei starting from re- 
alistic nucleon-nucleon (NN) and three-nucleon (NNN) 
interactions. On the other hand, a fully-developed fun- 
damental theory able to address a large range of nuclear 
scattering and nuclear reaction properties is still miss- 
ing, particularly for processes involving more than four 
nucleons overall. Better still, achieving a realistic ab ini- 
tio description of light nuclei requires abandoning the 
"traditional" separated treatment of discrete states and 
scattering continuum in favor of a unified treatment of 
structural and reaction properties. 

The development of such a unified fundamental theory 
is key to refining our understanding of the underlying 
forces across the nuclear landscape: from the well-bound 
nuclei to the exotic nuclei at the boundaries of stability 
that have become the focus of the next generation exper- 
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iments with rare-isotope beams, to the low-energy fusion 
reactions that represent the primary energy-generation 
mechanism in stars, and could potentially be used for 
future energy generation on earth. 

In the recent past, significant effort has been devoted 
to extend ab initio techniques to the treatment of dy- 
namical processes among light nuclei To this 
aim, we introduced a new many-body approach based 
on expansions over fully-antisymmetric {A — a, a) binary- 
cluster states in the spirit of the resonating-group method 
(RGM) [lol - [T5l |. in which each cluster of nucleons is de- 
scribed within the ab initio NCSM 16]. The unknown 
relative-motion wave functions between pairs of clus- 
ters are obtained by solving a set of non-local integral- 
differential coupled-channel equations and have appropri- 
ate bound-state and/or scattering asymptotic behavior. 
Capable of treating bound and scattering states of light 
nuclei in a unified formalism starting from the funda- 
mental inter- nucleon interactions, the NCSM/RGM ap- 
proach [13] has been successfully applied to a wide va- 
riety of binary processes, such as nucleon-^He and n— ^Li 
scattering [l^, ''Be(p,7)^B capture [l^, d— "^He scatter- 
ing d^, 3H(d,n)4He, and ^}ie{d,p)*lle fusion [U, and 
an extension to the treatment of three-cluster dynam- 
ics is under development [12, H^l- At the same time, 
these studies have highlighted practical limitations of the 
approach mainly related to a non-entirely efficient con- 
vergence behavior at short-to-medium distances, as dis- 
cussed in the following. 

Two kinds of convergence patterns have to be taken 
into account when performing a NCSM/RGM calcula- 
tion. First, one has to investigate the dependence on 
the size of the harmonic oscillator (HO) basis used to 
expand the NCSM eigenstates of the clusters and local- 
ized components of the couplings between binary-cluster 
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states. This size is characterized by A'max, the maxi- 
mal number of HO excitations above the lowest possi- 
ble configuration of the clusters. With soft similarity- 
renormalization- grou p (SRG) [23-127} evolved chiral EFT 
NN interactions [2l[29jl, employed in most NCSM/RGM 
calculations, HO basis sizes with iVmax ^ 10—14 are typi- 
cally sufficient to reach convergence and computationally 
feasible. Second, one has to study the convergence with 
respect to the number of clusters' eigenstates included in 
the calculation. While including only the ground state 
(g.s.) of the tightly-bound ^He in nucleon-^He scattering 
calculations already leads to a very good approximation 
of the A — 5 scattering phase shifts [l^, the descrip- 
tion of the low-energy ''Be(p,7)^B capture required tak- 
ing into account the lowest five eigenstates of ^Be 19]. 
The convergence with the number of clusters' eigenstates 
becomes even more problematic for weakly-bound clus- 
ters. Calculations with composite projectiles (the lighter 
of the two clusters) such as ^H, "^H, and '^He, show that 
it is essential to take into account the virtual breakup of 
these systems even at energies much below the breakup 
threshold. Presently, this is achieved by including a large 
number of excited pseudostates [13, [Hf, H^l of the projec- 
tile. This in turn results in a dramatic increase of com- 
plexity of the calculations as a large number of channels 
are coupled. 



In this paper we present a more efficient approach to 
nuclear bound and continuum states, the no-core shell 
model with continuum (NCSMC). We adopt an extended 
model space that, in addition to the continuous binary- 
cluster {A — a, a) NCSM/RGM states, encompasses also 
square-integrable NCSM eigenstates of the v4-nucleon 
system. Such eigenstates introduce in the trial wave 
function short- and medium-range A-nucleon correlations 
that in the NCSM/RGM formalism have to be treated 
by including a large number of excited states of the clus- 
ters. An analogous approach was suggested already in 
the original RGM papers [13, [HI ■ The idea behind the 
NCSMC was first mentioned in our review paper 13 ll 
and the formalism was succinctly introduced in Ref. [S^l , 
where it was applied to study of the low-lying resonances 
of the exotic ^He nucleus using an SRG-evolved chiral 
EFT NN potential that provides an accurate descrip- 
tion of the NN system. Here, we give a detailed presen- 
tation of the formalism, discuss the results published in 
Ref. 32], and present additional results. 



In Sec. ini we briefly review the NCSM and 
NCSM/RGM approaches and then introduce in detail 
the NCSMC formalism. In Sec. HIl we apply the NC- 
SMC to the exotic ''He nucleus. We discuss calculations 
presented in Ref. ISSj as well as additional results. Con- 
clusions and outlook are given in Sect. IIVI Parts of the 
formalism not suitable for the main text are presented in 
Appendix [X] 



II. FORMALISM 

This section is dedicated to the formalism of the NC- 
SMC theory with a particular focus on the case in which 
the binary-cluster portion of the basis is given by a single- 
nucleon projectile in relative motion with respect to an 
{A — l)-nucleon target. First, in Sec. Ill Al we briefly re- 
view the NCSM, then in Sec. Ill Bl we present useful back- 
ground and expressions for the NCSM/RGM formalism. 
Finally, in Sec. ITTCl we introduce in detafl the NCSMC. 

A. NCSM 

The ab initio NCSM is a structure technique appro- 
priate for the description of bound states or for approx- 
imations of narrow resonances. Nuclei are considered as 
systems of A non-relativistic point-like nucleons inter- 
acting through realistic inter-nucleon interactions, i.e., 
those that describe accurately two-nucleon and, possibly, 
three-nucleon systems. All nucleons are active degrees 
of freedom. Translational invariance as well as angular 
momentum and parity of the system under considera- 
tion are conserved. The many-body wave function is cast 
into an expansion over a complete set of antisymmetric 
A-nucleon HO basis states containing up to iVmax HO 
excitations above the lowest possible conflguration: 

N=0 i 

Here, N denotes the total number of HO excitations of 
all nucleons above the minimum conflguration, J'^T are 
the total angular momentum, parity and isospin, and i 
additional quantum numbers. The sum over N is re- 
stricted by parity to either an even or odd sequence. The 
basis is further characterized by the frequency of the 
HO well and may depend on either Jacobi relative or 
single-particle coordinates. In the former case, the wave 
function does not contain the center of mass (cm.) mo- 
tion, but antisymmetrization is complicated. In the lat- 
ter case, antisymmetrization is trivially achieved using 
Slater determinants, but the cm. degrees of freedom are 
included in the basis. The HO basis within the iVmax 
truncation is the only possible one that allows an exact 
factorization of the cm. motion for the eigenstates, even 
when working with single-particle coordinates and Slater 
determinants. Calculations performed with the two al- 
ternative coordinate choices are completely equivalent. 

Square-integrable energy eigenstates expanded over 
the A^niax^f^ basis, \ANiJ'^T), are obtained by diago- 
nalizing the intrinsic Hamiltonian, H = Ti„t + V, 

HIAXJ^'T) = Ex\A\rT) , (2) 

where Tint is the internal kinetic energy operator and V 
the NN or NN+NNN interaction. Convergence of the 
HO expansion with increasing A'max values is accelerated 
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by the use of effective interactions derived from the un- 
derlying potential model through either Lee-Suzuki simi- 



larity transformations in the NCSM space Il6l . 33 1 or SRG 
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^5 . In 



transformations in momentum space [241-127 
this latter case, the NCSM calculations are variational. 
Finally, we note that with the HO basis sizes typically 
used (A'lnax^lO— 14), the \AXJ'^T) eigenstates lack cor- 
rect asymptotic behavior for weakly-bound states and al- 
ways have incorrect asymptotic behavior for resonances. 



B. NCSM/RGM 

In the NCSM/RGM, the ansatz of Eq. (P) for the 
A-nucleon wave function is replaced by an expansion 
over antisymmetrized products of binary-cluster channel 
states l^';^^"^) ^^-d wave functions of their relative motion 



E 



dr r' 



(3) 



The channel states |<I>;^^-^) contain {A — a)- and a-nucleon 
clusters (with a < ^) of total angular momentum, parity, 
isospin and additional quantum number /i , tti , Ti , ai and 
^2, 7''2, T2, a2, respectively, and are characterized by the 
relative orbital angular momentum £ and channel spin 

s = Ii+l2: 
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X Yi{rA-a,a) 



(4) 



The channel index v collects the quantum numbers 
{A- a ail^^Ti^a a2l2^T2; si}. The intercluster relative 
vector rA-a,a is the displacement between the clusters' 
centers of mass and is given in terms of the single-particle 
coordinates fi by: 



1 



rA- 



^ A — a, a — a, a 
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A— a 
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i=l 



^ 1 
Ti 



A 

E ^ 

j=A-a+l 



(5) 

The cluster wave functions depend on translationally in- 
variant internal coordinates and are antisymmetric under 
exchange of internal nucleons, while the intercluster an- 
tisymmetrizer Ai, takes care of the exchange of nucleons 
belonging to different clusters. 

With appropriate boundary conditions imposed on the 
wave functions of the relative motion juir), the expan- 
sion of Eq. ([U is suitable for describing bound states, 
resonances and scattering states between clusters. For 
bound states, expansions ([1]) and ([3]) are equivalent, al- 
though for well-bound systems where short-range A-body 
correlations play a dominant role, the convergence of the 
eigenenergy would typically be more efficient within the 
NCSM model space defined by Eq. 

The unknown relative-motion wave functions (r) are 
determined by solving the many-body Schrodinger equa- 
tion in the Hilbert space spanned by the basis states 



-0, 



where 



(6) 

(7) 

(8) 



are the Hamiltonian and norm kernels, respectively, and 
E is the total energy in the cm. frame. 

When computing Eqs. ([7]) and ([5]), the "exchange" 
terms of the norm kernel arising from the non-identical 
permutations in Av as well as all localized parts of the 
Hamiltonian kernel are obtained by expanding the radial 
dependence of the basis states of Eq. Q on HO radial 
wave functions Rniir) according to: 
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Rneir)\<P 



(9) 



where P indicates the HO model space and 
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{\A - a aiIl'Ti)\a 02/^^2))^'^^ 



X Yi{fA-a,a) 



Rnt{rA-a,a) 



(10) 



Here, we remind that the A-nucleon microscopic 
Hamiltonian can be written in the form 



H = Trel + Kcl + Vc{r) + HiA-a) + ff(a) 



(11) 



where Tj-ei is the relative kinetic energy between target 
and projectile, Vrei includes all the interactions between 
nucleons belonging to different clusters after subtraction 
of the average Coulomb interaction between them Vc {r) 
(see [13 for a detailed discussion on this point), and 
H(A-a) and H(a) ^-re the intrinsic microscopic Hamilto- 



nians for A — a and a nucleons, respectively. The same 
inter-nucleon interactions are consistently employed in 
each term of Eq. pT|) . Accordingly, the clusters' eigen- 
states |A — a Qi/^^Ti) and \a Q!2/J^T2) are obtained 
by NCSM diagonalization of their respective microscopic 
Hamiltonians H(^A-a) and -ff(a)- The same frequency and 
consistent model-space size are used in the HO expan- 
sions of the clusters and localized parts of the integra- 
tion kernels. The size N^^x of the HO model space is the 
same for states of the same parity, whereas it differs by 
one unit for states of opposite parity. 

While the NCSM/RGM formalism has been fully de- 
veloped for single- (a = 1) [l3|, two- (a = 2) [S^l and 
three- nucleon (a = 3) projectiles [lOl, and can be also 
extended to a = 4 projectiles as well as to three-body 
clusters [2^ , in this work we limit ourselves to the a = 1 
case, where the inter-cluster antisymmetrizer is defined 
as 



(12) 
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and Pi^A is the permutation operator exchanging the i-th 
particle in the target with the projectile nucleon, labeled 
by the index A. 

1. Orthogonalizatton m the NCSM/RGM 

Here, we recall some of the details concerning the or- 
thogonalization of the NCSM/RGM equations © that 
are useful for our further discussion of the NCSMC for- 
malism. 

Because of the non-identical permutations in the inter- 
cluster antisymmetrizer, the channel states Av\'^il'^) are 
not orthonormal to each other. In general, we prefer to 



work with the orthonormalized binary-cluster states 

J dr'r'' N;) (r, r') AA<:J) , (13) 

where we introduced the inverse square root of the 
NCSM/RGM norm kernel In the following we re- 
view how this as well as the square root of the norm 
kernel are obtained. 

As anticipated in the previous section, the "exchange" 
term arising from the permutations in that differ 
from the identity are obtained using the HO expansion 
of Eq. Hence, using Eqs. ([T^ and ([TU|. the r-space 
representation of the norm kernel can be written as 



AC7(r, r') = 5,J-^^^^ - - 1) E Rni{r){'^i7\pA-,A^iZ)Rn'e{r') 

n,n' 

nn' 



6{r-r') 



nn' GP 



where we introduced the model-space norm kernel: 

^u2'n- = ^i^u'Snn' ^ - 1) ^ Rnl{r) Rn' t' {r' ) 

n.7i'GP 

X {<i>rj\pA-uA\^i:,^) ■ (15) 



The last line of Eq. (|T4|) shows that the r-space repre- 
sentation of the kernel is given by the convolution of the 
model-space kernel (second term) plus a correction due 
to the finite size of the model space P (first term). Square 

and inverse-square roots Af^^^f (r, r') can then be defined 
in an analogous way as: 



= 5vw' 



5{r — r') ■r-^ , 

2^ Rndr)dnn'Rn't{r ) 



nn' 

±1 



(16) 



is obtained from the hermitized Hamiltonian kernel, 



A-a A-a 



(18) 



i i 

for which we have borrowed the same notation H^^i (r, r') 
used previously in Eq. ([?])• 

Finally, the orthogonalized RGM equations read 



(19) 



with the wave functions of the relative motion Xiy(r) re- 
lated to the original functions ^^{r) by 



7/' 



(20) 



nn' 



For more details on the NCSM/RGM kernels we refer the 
interested reader to Ref. [17i] . 



where the model-space square and inverse square roots 
^vnv'n' obtained from the spectral theorem. 

The NCSM/RGM Hamihonian kernel within the or- 
thonormal basis of Eq. 



H..'(r, r') 



(17) 



C. NCSMC 

The NCSMC ansatz for the many-body wave function 
includes both A-body square-integrable and {A — a, a) 
binary-cluster continuous basis states according to: 



X 



|AAJ-T)+E / drr^^-^AA<^) 



(21) 
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The resulting wave function (|2T|) is capable of describing 
efficiently both bound and unbound states. Indeed, the 
NCSM sector of the basis (eigenstates \ AXJ'^T)) provides 
an effective description of the short- to medium-range 
A-body structure, while the NCSM/RGM cluster states 
make the theory able to handle the scattering physics of 
the system. In other words, with the expansion ([21]) one 
obtains the coupling of the NCSM with the continuum. 
Clearly, the NCSMC model space is overcomplete, but 
this is not a concern, as it will be shown in the following. 



1. NCSMC equations 

The discrete (ca) and continuous (7i/('")) unknowns of 
the NCSMC wave function are obtained as solutions of 
the coupled equations 



Hncsm h_ 

h n 



= E 



1 9 
9 1 



(22) 



where Xnir) Eire the relative wave functions in the 
NCSM/RGM sector when working with the orthogonal- 
ized cluster channel states of Eq. p^ . These are related 
to the original wave functions ^i, (r) of Eq. (PT|) by the re- 
lationship (l20t . Note, however, that the Xv{f) appearing 
in Eqs. and (|22l) are in general different, i.e., they 
are solutions of different equations. 

The NCSM sector of the Hamiltonian kernel is a diag- 
onal matrix of the NCSM energy eigenvalues E\ ([2]), 



NCSM)\X 



)aa' = {A\J^T\H\A\'J^T) = E),S, 



(23) 



while H is the orthogonalized NCSM/RGM kernel of 
Eq. (fT7|) . Because of the orthogonalization procedure of 
Sec. IIIB 1[ both diagonal blocks in the NCSMC norm 
kernel N are identities in their respective spaces 



N, 



XX' 



Sxx' gxu'{r') 
9X'u{r) 5uu'^^~p^ 



(24) 



The coupling between square-integrable and binary- 
cluster sectors of the model space is described by the 
cluster form factor 



in the norm kernel, and by the coupling form factor 



(25) 



hxu{r) 



dr'r'\A\rT\HAAK'r' ) K~i (^', 



(26) 

in the Hamiltonian kernel. Detailed expressions for these 
form factors are given in Appendix [A) The calculation 
of {A\J'^T\Au'^'l'l'^) overlap matrix elements between 
NCSM wave functions and binary-cluster states was also 
discussed in Ref . [s^ . We also note that by squaring the 
absolute value of these matrix elements and integrating 
over r, one obtains spectroscopic factors. 



The NCSMC equations can be orthogonalized in an 
analogous way to that presented for the NCSM/RGM in 
Sec. Ill B II To define the square and inverse square root 
of the NCSMC norm in the r-space representation, we 
first rewrite Eq. (|24| as the convolution of the model- 
space norm kernel plus a correction for the finite size of 
the HO model-space P 



N, 



XX' 



Sxx 

Ri^rC 



N, 



XX' 



'X'X' " 

Ry' r' u' 7l' 



(27) 



where the sum over the repeating indexes X,i',n, X' , D' , 
and n' is implied, the notation RurOn stands for 
Rne{r)6^o, and the model-space NCSMC norm is given 
by: 



Sxx 



XX' dxO'n' 
9x'i>n Si^O'Snn' 



(28) 



Here, the model-space cluster form factor is related to the 
r-space one through gx^ir) = J2n Rni{r)gxiyn (as demon- 
strated in Appendix]^. Accordingly, the square and in- 
verse square roots of N can then be defined as: 

fiV±l)AA', , 





5uu' ^'-''SlJ •* - Rni{r)5u^'5nn'Rn't{r') 



^XX ] /^±i^AA• / -A'A 
R^ru ' ^ '-^n^'n' \ n 



h'X' 
Ru' r' D' n' 



(29) 



Inserting the identity A^^^TV+I in both left- and right- 
hand sides of Eq. ((22|) and multiplying by N^^s from the 
left one obtains the orthogonalized NCSMC equations 



where the orthogonalized Hamiltonian is given by. 



(30) 



and the orthogonal wave functions by: 



X 



X 



(32) 



Finally, the ansatz (l^lT) in terms of the orthogonalized 
NCSMC wave function takes the form: 
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+ J2 drr' dr' r'^AA^f/W;) (r, r') 



.(33) 



2. Solving the NCSMC equations 

At large inter-cluster distances r, the clusters are as- 
sumed to interact through the Coulomb interaction only. 
Hence, the NCSMC equations are solved dividing the 
space into an internal region, r ^ tq, and an external 
region, r > rg, and applying the coupled-channel micro- 
scopic R-matrix method on a Lagrange mesh The 
separation radius r — must be large enough to ensure 
that the wave function of the A-body states \AXJ'^T) 
vanishes when approaching the external region, where 
the asymptotic behavior of the NCSMC solutions is de- 
scribed entirely by the radial wave functions 



(34) 



and 



(35) 



for bound and scattering states, respectively. Here, 
Wi{r]^, K^r) are Whittaker functions and Hj^{ri^, k^t) are 
the incoming and outgoing Coulomb functions, with 
the speed, k^, the wave number, and 77^ the Sommerfeld 
parameter of the final state being studied. Asymptotic 
normalization constant for bound states and scattering 
matrix between initial (j) and final {v) scattering states 
are denoted respectively with ^ and S^l^ . The func- 
tions wl ^(r) stand for either the non-orthogonalized 
wave function Xi^(j') or for the orthogonalized Xi^{r) ac- 
cording to which set of equations, Eq. (1^ or ([50]) . is 
being considered. 

One of the advantages of the microscopic R-matrix 
method is that the wave function in the internal re- 
gion can be expanded on a set of square-integrable func- 
tions. A particularly convenient choice when dealing 
with non-local potentials, as in our case, is the set of 
Lagrange functions /n(r) associated with the shifted 
Legendre polynomials and defined on the mesh points 
r„ e (0, To) [33|, labeled by the index 1 < n < iV. When 
the Gauss quadrature approximation is adopted, the La- 
grange functions are orthogonal to each other. In addi- 
tion, thanks to the Gauss quadrature approximation and 
the properties of the Lagrange functions, matrix elements 
of non-local potentials are proportional to the values of 
the non-local potentials at the mesh points. The number 
of mesh points N_ has to be large enough to guarantee 



an accurate representation of the wave functions in the 
internal region up to the matching radius tq. Typically, 
25 mesh points are sufficient to calculate a phase shift 
within six significant digits for tq — 15 fm. 

The matching between internal and external regions, 
and hence the imposition of the asymptotic behavior of 
Eqs. (p4)) and/or ([35| . is ensured by the Bloch surface 
operator (here generalized to account for the A-body 
square-integrable sector of our basis) 







A^5(r-a)(A_S.) 



and solving the Bloch-Schrodinger equations 



{H + L-E) 



X 



(36) 



(37) 



The operator H + L is Hermitian when the boundary 
parameter B^, is real. Because of the Bloch operator, 
the wave function in the right-hand side of Eq. (j37p can 
be replaced by its asymptotic behavior. When searching 
for bound states, B^, is chosen in such a way that the 
right-hand side vanishes, and one is left with the diago- 
nalization problem: 



(H + L) 



= E 



(38) 



For scattering states, the scattering matrix and the scat- 
tering wave functions are computed by solving Eq. ([57)1 
with the boundary parameter B^ — for each value of 
the relative kinetic energy E^in of the projectile-target 
system. The phase shifts S{Ekin) can then be extracted 
from the S-matrix. Energetically open and closed chan- 
nels are treated on equal footing. 



III. APPLICATION TO ^He 

The ^He nucleus is a particle-unstable system with a 
J'^T = 3/2- 3/2 ground state lying at 0.430(3) MeV 
[Hillil above the ^He-|-n threshold and an excited 5/2~ 
resonance centered at 3.35 MeV, which mainly decays to 
a + 3n (as discovered in the pioneering work of Ref. pO|). 
While there is a general consensus on the 5/2~ state, dis- 
cussions are still open for the other excited states. In par- 
ticular, the existence of a low-lying {En^l MeV) narrow 
(r < 1 MeV) l/2~ state has been advocated by many 
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experiments [4l|-|45j (most of them using knockout reac- 
tions with a ®He beam on a carbon target), while it was 
not confirmed in several others [46l - [5ll |. This contradic- 
tory situation arises from the main experimental diffi- 
culty of measuring the properties of ^He excited states 
in the presence of a three-body background of ^Hc plus 
n (coming from the particle decay of ''He) plus a third 
outgoing particle involved in the reaction used to pro- 
duce ^He. In addition, as it has been pointed out in 
one of the most recent experimental works [Slj , some of 
the earlier data could have been affected by background 
noise coming from the interaction with the carbon tar- 
get, while a polypropylene (CH2)n target would reduce 
the background contamination. The presence of a low- 
lying l/2~ state has also been excluded at the 90% con- 
fidence level by a study on the isobaric analog states of 
^He in ^Li 52|. According to this work, a broad 1/2^ 
resonance at ~3.5MeV with a width F'-^lOMeV fits the 
data the best. Neutron pick-up and proton- removal reac- 
tions [11,14^ suggest instead a l/2~ resonance at about 
3 MeV with a width T « 2 MeV. 

From a theoretical point of view, ''He is an ideal system 
to showcase new achievements made possible by a uni- 
fied ah initio approach to nuclear bound and continuum 
states such as the NCSMC. Since ^He is unbound, it can- 
not be reasonably described within the NCSM. One could 
calculate its properties using the NCSM/RGM within an 
^He-|-n binary-cluster expansion. However, the ^He nu- 
cleus is weakly bound and all its excited states are un- 
bound. Consequently, a limitation to just a few lowest 
^He eigenstates in the NCSM/RGM expansion would be 
questionable especially because, except for the lowest 2+ 
state, all other ^He excited states are either broad res- 
onances or simply states in the continuum. As we will 
show in the following, with the NCSMC these challenges 
are overcome. Finally, for this study we use the SRC 
evolved [2i-'27l chiral N3L0 NN potential (500 MeV cut- 
off) of Refs. (28fci29||. For the time being, the induced and 
initial chiral three-nucleon interactions arc not included 
in the calculations, therefore our results depend on the 
low-momentum SRC evolution parameter A. However, 
by selecting A — 2.02 fm^^, we obtain very realistic 
binding energies for the lightest nuclei, e.g., the ^He (see 
Table [l| and, more importantly for the present investi- 
gation, the ^He. Consequently, this choice of NN poten- 
tial allows us to perform qualitatively and quantitatively 
meaningful calculations for ''He that can be compared 
to experiment. In the following sections, we discuss the 
convergence of the NCSMC calculation and compare it 
to the corresponding NCSM and NCSM/RGM results. 
We will also address the controversial issue of a low-lying 
l/2~ resonance in ^He. 



E^.s. [MeV] 


*He 


"He 


7tt 

^He 


NCSM iVn.ax = 12 


-28.05 


-28.63 


-27.33 


NCSM extrap. 


-28.22(1) 


-29.25(15) 


-28.27(25) 


Expt. 


-28.30 


-29.27 


-28.84 



TABLE I: Ground-state energies of "'^•^He in MeV. NCSM 
calculations were performed with the SRG-N'^LO A^A'' poten- 
tial with A = 2.02 fm-^ The HO frequency m=l& MeV 
was used in the shown A^niax=12 calculations. Exponential 
extrapolation was employed. 



input for the subsequent NCSAI/RGAI and NCSMC in- 
vestigations of ''He. 

Our calculated ^He ground-state energies for a range 
of HO frequencies and various basis sizes (A^max values) 
are presented in Fig. [T] The variational NCSM calcula- 
tions converge rapidly and can be easily extrapolated to 
A^max oo using, e.g., an exponential function of the 
type -E(A^max) — E^o + a e"^^'"'"'. Results of such an ex- 
trapolation are shown in Fig.[2]where the A^max=8, 10 and 
12 points were used to determine the fitting parameters. 
The extrapolated ground-state energy with its error es- 
timate, based on extrapolations at other frequencies and 
different point selections, and the calculated energy at 
N^^^ = 12, nn = 16 MeV are given in Table HI The 
calculated ^He g.s. energy agrees quite well with experi- 
ment on both the absolute value and the separation with 
respect to the ^He-|-2n threshold. 

As shown in Fig.[Tl at A^^ax = 12 the dependence of the 
^He g.s. energy on the HO frequency is flat in the range 
of HQ, ~ 16 — 19 MeV. In general, when working within 
an HO basis, lower frequencies are better suited for the 
description of unbound systems. Therefore, we choose 



-12 
-14 
-16 
-18 
> -20 
S -22 
«^"-24 
-26 
-28 
-30 



1 — ' — I — ' — I — ' — r 



T — ' — I — ' — r 

He SRG-nYo NN 



A=2.02 fm 




















A — AN^„=10 






--V - 


+ . . + extrap 





— ^ — I — ^ 



J 1 I 1 L 



J 1 I 1 I 1 L_ 



12 13 



14 



15 16 17 18 19 
hSi [MeV] 



20 



21 22 



A. "^He and ^He NCSM calculations 

We begin our discussion of results with NCSM calcu- 
lations for ^He that will generate eigenstates needed as 



FIG. 1: (color online). Ground-state energy of ®He calculated 
within the NCSM using the SRG-N^LO NN potential with 
A = 2.02 fm~^. The dependence on the HO frequency for 
different Amax basis sizes is shown. The points with error 
bars represent the results of the exponential extrapolation. 
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FIG. 2: (color online). Basis size Nmax dependence of ®He 
and ^He ground-state energies calculated within the NCSM. 
The SRG-N^LO NN potential with A = 2.02 frn'^ and the 
HO frequency of hQ,=lQ MeV were used. Exponential extrap- 
olations from the last three A^max points is shown. 



hSl = 16 MeV for the calculation of the other ^He eigen- 
states that will be used as input for the NCSM/RGM 
and NCSMC investigations of the ^He nucleus. At the 
same time, we also performed NCSM/RGM and NCSMC 
calculations with h£l — 19 MeV, to test the stability of 
our results against this parameter. Calculated ^He exci- 
tation energies for basis sizes up to A'max = 12 are shown 
in Fig. [H The ^He nucleus is a is a weakly-bound Bor- 
romean system. All its excited states are unbound, and, 
except for the lowest 2+, either broad resonances or states 
in the continuum. The excitation energy of the 2^ state 
is fairly stable with respect to the basis size of our NCSM 
calculations. The higher excited states, however, drop in 
energy with increasing N^^y^ with the most dramatic ex- 
ample being the multi-particle-hole Oj^ state. This spells 
a potential difficulty for the NCSM/RGM calculations 
as, with increasing density of ^He states at low energies, 
a truncation of the model space to include just the few 
lowest cigcnstates becomes questionable. In addition, in 
a NCSM/RGM study of '^He one should also consider the 
contribution of binary-cluster states in which the neutron 
is coupled to negative-parity states of ^He, where similar 
issues arise. 

Next we performed NCSM calculations for ^He ground 
and excited states, which will serve as input to the NC- 
SMC calculations described in the next section. The cal- 
culated g.s. energy for different basis sizes is shown in 
Fig. [5] together with the exponential extrapolation and 
the ^He g.s. energies discussed earlier. The largest-space 
values and the extrapolated energies are also given in 
Table m The NCSM calculation predicts ^He unbound in 
agreement with experiment. However, the resonance en- 
ergy with respect to the ^He+n threshold appears overes- 
timated (contrary to the ^He^"'lle+2n case). Obviously, 
it is not clear that the ad hoc exponential extrapolation 
is valid for the unbound states. Or, it may have a siz- 
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FIG. 3: (color online). Dependence of ^He excitation ener- 
gies on the size of the basis Nmax- NCSM calculations were 
performed with the SRG-N^LO NN potential with A = 2.02 
fm~^ and the HO frequency of hfl=16 MeV. 



able systematic uncertainty compared to the bound-state 
case. Nevertheless, the differences between iVmax = 12 
and extrapolated energies suggest that the fastest con- 
vergence rate is obtained for the strongly-bound "^He and 
the slowest for the unbound ''He, as one would expect. 
Finally, we note that no information on the width of the 
resonance can be obtained from the NCSM calculation, 
which is performed in a square-integrable HO basis. 

However, we can study the structure of the ^He NCSM 
eigenstates by evaluating their overlap functions with 
^He-|-n binary-cluster channels. These overlap functions, 
or cluster form factors, g\^{r) [see Eqs. ([23, ((X3|) ] 
are also one of the inputs to the NCSMC calculations. 
By integrating g^jyir) over r, we obtain the spectroscopic 
factors summarized in Table HIl Note that there we use an 
alternative coupling scheme [compared to Eq. ^] more 
commonly used in the literature for spectroscopic fac- 
tors. Overall, we find a very good agreement with the 
variational Monte Carlo (VMC) and GFMC results as 
well as with the latest experimental value for the ground 
state [111 . Interesting features to notice are the spread of 
the 3/2~ g.s. wave function over all three considered ^He 
states with a dominance of the 2f and the about equal 
spread of l/2~ ^He excited state between the ^He 0+ and 
2^ states. The ^He 5/2~ state has about the same con- 
tributions from the 2f and 2j ^He states with the former 
of an almost pure s=5/2 component (with s the channel 
spin defined in Eq. Though spectroscopic factors 

are not observable, they provide valuable information on 
the structure of the wave function. In the present study, 
overlap functions and spectroscopic factors are not the 
final products to be compared to experiment, but rather 
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3/2r 


2+ 




0.42 




0.30 




1/2- 


0+ 




0.94 


0.69 


0.91 
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0.34 
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0.26 
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0.93 
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0.81 




5/2- 


2t 


-Pi 


0.49 
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0.37 
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2^ 


-Pi 


1.10 


1.05 


1.07 




3/22 
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-Pi 


0.08 


0.32 


0.03 




3/22- 


2+ 

^2 


-p\ 


0.03 








3/22- 


2+ 

^2 


-Pi 


0.25 









TABLE II: NCSM spectroscopic factors com par ed to Cohen- 
Kurath (CK) JSl] and VMC/GFMC [J, [H,!!! calculations 
and experiment. NCSM calculations were performed with the 
SRG-N=*LO potential with A=2.02 fm-\ N^^^=V1 and 
the HO frequency of ^^=16 MeV. The CK results should in 
principle be still multiplied by j4/(yl— 1) to correct for the 
center of mass motion. 



inputs to more sophisticated NCSMC calculations. 
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FIG. 4: (color online). Dependence of the NCSM/RGM (a) 
and NCSMC (b) ^He + n phase shifts of the ''He 3/2- ground 
state on the number of ®He states included in the binary- 
cluster basis. The short-dashed green curve, the dashed blue 
curve and the solid red curve correspond to calculations with 
®He 0"*" ground state only, 0^, 2"*" states and 0"*", 2^, 2"*" states, 
respectively. The SRG-N^LO iViV potential with A = 2.02 
fm^'^, the A'^max=12 basis size and the HO frequency of ftf2=16 
MeV were used. See text for further details. 



B. 'He NCSM/RGM and NCSMC calculations 

In the following, we present NCSMC calcula- 
tion for the ^He nucleus performed within a model 
space containing the six lowest negative-parity 
(3/2];, 1/2-, 5/2-, 3/2^, 3/2^, 3/24) and four lowest 
positive-parity (1/2+, 5/2+, 3/2+, 5/2+) NCSM eigen- 
states of ''He plus n+^He NCSM/RGM binary-cluster 
channels including up to the three lowest eigenstates of 
^He, i.e.. 0+,2+, and 2^. For the sake of comparison, 
we will also present results obtained by retaining only 
the binary-cluster portion of such a model space [i.e., 
only the second term in Eq. (PT|) or, equivalently, the 
ansatz ^\ and solving the orthogonahzed NCSM/RGM 
equations of Eq. (fT9|) . 

We start by studying the dependence of the 3/2" g.s. 
phase shifts on the number of ^He eigenstates included in 
the NCSM/RGM [panel (a)] and NCSMC [panel (b)] cal- 
culations, shown in Fig. ID Here, the channels are denoted 
using the standard notation ^''+^£j, e.i?., ^^3/2 for the g.s. 
resonance, with the quantum numbers s, i and J defined 
as in SecHTH Eq. We observe that the NCSM/RGM 



calculation with the ^He target restricted to its ground 
state does not produce a ''He 3/2- resonance (the phase 
shift does not reach 90 degrees and is less than 70 de- 
grees up to 5 MeV). A ^^3/2 resonance does appear once 
n-|-^He(2+) channel states are coupled to the basis, and 
the resonance position further moves to lower energy with 
the inclusion of the second 2+ state of ^He. On the con- 
trary, the NCSMC calculation with only the ground state 
of ^He already produces the ^^3/2 resonance. In fact, this 
NCSMC model space is sufficient to obtain the ^He 3/2- 
g.s. resonance at about 1 MeV above threshold, which 
is lower than the NCSM/RGM prediction of 1.39 MeV 
when three ^He states are included. Adding the first 2+ 
state of ^He generates a modest shift of the resonance 
to a still lower energy while the 2^ state of ^He has no 
significant infiuence [see Fig. 21 panel (b)]. We further 
observe that the difference of about 0.7 MeV between 
the NCSM/RGM and NCSMC results for the resonance 
position is due to additional correlations in the wave func- 
tion brought about by the 'He eigenstates that are cou- 
pled to the neutron-^He binary-cluster states in the NC- 
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FIG. 5: (color online). Dependence of the NCSM/RGM (a) 
and NCSMC (b) ^He + n phase shifts of the '^He 1/2" excited 
state on the number of ®He states included in the binary- 
cluster basis. The short-dashed green curve, the dashed blue 
curve and the solid red curve correspond to calculations with 
^He 0^ ground state only, 0^, 2"*" states and 0"*", 2^, 2^ states, 
respectively. See Fig. U for further details. 

SMC. Indeed, such A — 7 eigenstates (in the present 
calculation four 3/2" states, of which only the 3/2^ pro- 
duces a substantial effect on the ^^3/2 resonance) have 
the practical effect of compensating for higher excited 
states of the ^He target omitted in the NCSM/RGM 
sector of the basis. These omitted ^He states include 
both positive-parity, some of which are shown in Fig. [31 
and negative-parity excitations such as, e.g., the 1~ soft 
dipole excitation etc. While NCSM/RGM calculations 
with a large number of excited states of the target or 
projectile can become prohibitively expensive, the cou- 
pling of a few square-integrable NCSM eigenstates of the 
composite system is straightforward. Because of this, the 
NCSMC approach offers a superior rate of convergence 
and is much more efhcient, as demonstrated in Fig. 31 

A similar, although less dramatic, difference between 
NCSM/RGM and NCSMC calculations is shown in Fig. El 
for the 1/2" excited state of ^He. Here, the res- 
onance is quite broad with a slowly increasing phase 
shift. It is interesting to note that the 1/2" state couples 
strongly to the 2^ state of ^He (the spectroscopic factor 
is large, see Table IIIl) . This causes a small but visible shift 
of the phase shift when this state is added to the 



FIG. 6: (color online). Dependence of the NCSMC ''He + n 
phase shifts of ^He 3/2" (a), 1/2" (a) and 5/2" (b) states on 
the size of the HO expansion Nmax- The ^He 0"*", 2+, 2+ states 
were included in the binary-cluster basis. The SRG-N'^LO 
A^A'^ potential with A = 2.02 fm"^ and the HO frequency of 
hQ—16 MeV were used. 



NCSMC calculation in panel (b) (full vs. dashed line). 
The 1/2" state presents a significant overlap also with 
the 0+ and 2f states of ^He, and this is the reason of its 
broadness. 

In Fig. |6l we present the dependence of the NCSMC 
^P3/2, [panel (a)] and ^^5/2 [panel (b)] phase shifts 

on the size of the HO basis in the range 6 < A^^max < 
12. While the variation between A^max = 6 and Amax = 
8 curves is substantial, it becomes quite small between 
A^max = 10 and A^max = 12 results. Based on this, we do 
not expect that an A^max = 14 calculation, which at this 
time is computationally out of reach, would significantly 
change the present A^ax — 12 picture. 

The NCSM/RGM and NCSMC phase shifts for the 
71 + ^He five P-wave and the ^5*1/2 channels are shown in 
Fig. [71 All curves have been obtained including the three 
lowest ^He states (i.e., the 0+ ground state and the two 
lowest 2+ excited states) within the A'niax = 12 HO basis. 
The model space of the NCSMC calculations [panel (b)] 
additionally includes ten ^He NCSM eigenstates, as de- 
scribed at the beginning of this section. The dashed verti- 
cal area centered at 0.43 MeV indicates the experimental 
centroid and width of the ^He ground state (38l. [39j. As 
expected from a variational calculation, the introduction 
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FIG. 7: (color online). NCSM/RGM (a) and NCSMC (b) 
®He + n diagonal phase shifts (except ^^3/21 which are eigen- 
phase shifts) as a function of the kinetic energy in the center of 
mass. The dashed vertical area centered at 0.43 MeV indicates 
the experimental centroid and width of the ^He ground state 
[3^. [33]. In all calculations the lowest three ^He states have 
been included in the binary-cluster basis. The SRG-N^LO 
NN potential with A = 2.02 fm"^ wfthin the iVmax = 12 ba- 
sis size and the HO frequency of HQ— 16 MeV were used. See 
text for further details. 



of the additional square-integrable A-body basis states 
I^AJ'^r) [i.e., going from panel (a) to panel (b) of Fig. [7] 
lowers the centroid values of all ^He resonances. In par- 
ticular, the 3/2~ ground and 5/2~ excited states of ''He 
are pushed toward the ^He -I- n threshold, closer to their 
respective experimental positions. The resonance widths 
also shrink toward the observed data as we discuss below. 
We note that we also calculated higher partial waves, e.g. 
ZJ-waves, in both approaches. However, the correspond- 
ing phase shifts are very small and do not present any 
interesting structures in the energy range displayed in 
Fig. [71 Therefore, we did not include them in the figure. 
Unlike the P-wave resonances, the influence of the ^He 
positive-parity NCSM eigenstates on these phase shifts 
is rather weak. 

The experimental centroid of the accepted 3/2^ and 
5/2~ resonances in ^He as well as the possible l/2~ state 
at 3.03 MeV [i^ are shown in Fig. [Altogether with our 
Nma,^ = 12 predictions. For NCSM/RGM and NCSMC, 
the resonance centroids are calculated as the values of 
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FIG. 8: (color online). Experimental and theoretical centroid 
energies for ■^He resonances, with the ^He -I- n threshold as 
the energy reference. The experimental energy of the l/2~ 
resonance is taken from Ref. [ig]. The theoretical values for 
NCSM/RGM and NCSMC correspond to the n-'Re kinetic 
energy in the center of mass when the derivative of the phase 
shift is maximal, see text for details. The information on the 
width of the states is given in Table IIIII The calculations are 
carried out as described in Tab.[Tl Fig. [7] and in the text. 



the kinetic energy in the center of mass Ekin for which 
the first derivative of the phase shifts is maximal [STj . 
The resonance widths are subsequently computed from 
the phase shifts according to (see, e.g., Ref. |58|): 



r = 



dS{Ek^n)/dE 

kv 



(39) 



Ekin — Ei^ 



where Er is the resonance centroid, evaluated as dis- 
cussed above, and the phase shift are expressed in radi- 
ans. Computed widths and Er values are reported in 
Table IIIII together with the available experimental data. 
An alternative, though less general, choice for the reso- 
nance energy could be the kinetic energy corresponding 
to a phase shift of 7r/2 (dashed horizontal lines in Fig. [7]). 
While the procedure of Eq. is safely applicable to 
sharp resonances, broad resonances would in principle re- 
quire an analysis of the scattering matrix in the complex 
plane. Here, we are more interested in a qualitative dis- 
cussion of the results, and will use the above extraction 
procedure for broad resonances as well. Though the two 
alternative ways of choosing Er lead to basically identi- 
cal results for our calculated 3/2 resonance, the same 
is not true for the broader 5/2^ resonance and the very 
broad l/2~ resonance. The less general 7r/2 condition, 
which is not valid for broad resonances, would result in 
Eji-^?,.7 MeV and r - 2.4 MeV for the 5/2" resonance 
and ERr^ A MeV (see Fig. [7]) and F - 13 MeV for the 
l/2~ resonance. 



12 



r 


experiment 




NCSMC 


NCSM/RGM 


NCSM 




Er V 
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Er V 
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Er 


3/2- 


0.430(3) 0.182(5) 


[39] 


0.71 0.30 


1.39 0.46 
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5/2' 
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[59] 


3.13 1.07 


4.00 1.75 
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3.03(10) 2 


m] 


2.39 2.89 


2.66 3.02 


3.26 




3.53 10 


[52] 










1.0(1) 0.75(8) 


[42] 









TABLE III: Experimental and theoretical values for the reso- 
nance centroids and widths in MeV for the 3/2' ground state 
and the 5/2- and 1/2- excited states of ^He. Calculations 
are carried out as described in Tab. [H Fig. [7| and in the text. 




Interestingly, the NCSM eigenenergy for the 3/2]^ 
ground state resonance is close to the energy centroid 
found within the NCSM/RGM approach. This is ac- 
cidental as both calculations are deficient in different 
ways. The NCSM lacks the description of long-range 
correlations due to the HO basis truncation, while the 
NCSM/RGM lacks a proper description of short- and 
medium-range correlations due to the omission of higher 
excited states of the ^He target. In the NCSMC, a sig- 
nificant energy shift is brought by the coupling of the 
two basis, with a quenching of the separation energy by 
almost 0.7 MeV, closer to the experimental findings. At 
the same time, in the NCSMC calculation, the resonances 
become sharper, with narrower widths, once again in a 
better agreement with experiment. Our NCSMC 3/2^ 
g.s. resonance position and width slightly overestimate 
measurement (e.g., the latest determination from the re- 
coil proton tagged knockout reaction for ^He [s^ finds 
Er = 0.430(3) MeV and T = 0.182(5) MeV). At the 
same time, predictions for the 5/2" resonance are lower 
compared to experiment , though our determina- 

tion of the width should be take with some caution in 
this case. 

In all three approaches considered here, the 1/2' res- 
onance is predicted below the 5/2" excited state. At the 
same time one has to keep in mind that the NCSM ap- 
proach is not expected to provide a reliable description of 
broad resonances and that our determination of the 1/2' 
resonance position in the NCSM/RGM and NCSMC has 
to be taken with some caution, as explained above. From 
an experimental standpoint, the situation concerning the 
1/2" resonance is not clear as discussed in the begin- 
ning of this section and documented in Table [Till While 
the centroid energies determined in the experiments of 
Refs. [H,!!^ and (s^l are comparable, the widths are very 
different. Within the present determination of Er and F, 
the NCSMC results are in fair agreement with the 1/2' 
properties measured in the neutron pick-up and proton- 
removal reactions experiments of Refs. [48| and [43|. Our 
calculations definitely do not support the hypothesis of 
a low-lying (^i^~l MeV) narrow (F < 1 MeV) 1/2' res- 
onance l4l|-[45l. 



FIG. 9: (color onhne). NCSMC ''He-hn 3/2- P-wave eigen- 
phase shifts as a function of the kinetic energy in the center 
of mass. Calculations are carried out as described in Fig. [7] 
(b). See text for further details. 

We also note that our NCSMC calculations predict two 
broad ^^3/2 resonances (dominated, respectively, by the 
first and second 2+ states of ^He) at about 3.7 MeV and 
6.5 MeV with widths of 2.8 MeV and 4.3 MeV, respec- 
tively. As shown in Figs. [7] and [HI the corresponding 
eigenphase shifts do not cross 7r/2. In Fig.[9l we present 
all P-wave eigenphase shifts in a broader energy range 
up to 10 MeV. There is a considerable mixing of the P 
waves around the 3/2^ resonance as it can be seen by 
comparing the eigenphase shifts of Fig. [9] with the diag- 
onal and phase shifts of Fig. [TUl The mixing 
parameter for other resonances is very small (of course, 
there is no mixing below the n-|-^He(2^) state threshold). 
In experiment, there is a resonance of undetermined spin 
and parity at 6.2(3) MeV with a width of 4(1) MeV [soj . 

The level order predicted in other theoretical calcula- 
tions mostly agrees with our present findings (T], [g^] . 
The widths of the ^He states were calculated recently in 
a "'He+n + n + n cluster model jg^l- The 1/2- state 
was found at low excitation energy (^1.05 MeV), but 
with a width of 2.19 MeV, i.e., close to what we find. 
The width of the 5/2" resonance, 1.5 MeV, obtained in 
Ref. [S^l is also comparable to our prediction. Two 3/2- 
resonances in addition to the g.s. resonance were reported 
in Ref. [6^. One of them just above the 5/2- state with 
a width of 1.95 MeV, while the other at the excitation en- 
ergy of about 5.3 MeV and a width of 5.77 MeV. This is 
qualitatively similar to our results although, in our case, 
the 3/22^ resonance is broader by 0.85 MeV. 

The ^He resonances were also investigated in RGM cal- 
culations of Ref. [6l| using a semi-realistic NN potential. 
The ordering of the resonances found in this study is the 
same as in ours and the ^^3/2 g-S. resonance phase shift 
is also in close agreement with our results. On the other 
hand, unlike in our calculations, the 1/2- and 5/2- P- 
wave resonances of Ref. [61!] do not cross 7r/2. While the 
^^3/2 and ^^3/2 diagonal phase shifts qualitatively agree 
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FIG. 10: (color online). NCSMC "^He + n 3/22 and 5/2" 
diagonal P-wave phase shifts as a function of the kinetic en- 
ergy in the center of mass. The calculation is carried out as 
described in Fig. [7](b). See text for further details. 



with ours, interestingly, the 5/2 resonance appears in 
the ^^5/2 partial wave rather than in the ^P^|2^ as found 
in our_calculations. The 5/2 ~ P- waves are reversed in 
Ref. 1^ compared to our calculations (see Fig. [TU] for 
the ^^3/2, ^Pzji, ^P^i2^ and ^^5/2 diagonal phase shifts). 

The Helium isotope g.s. properties, including those 
of ''He, were also recently investigated within the com- 
plex coupled-cluster method [gll- Using a realistic low- 
momentum NlSi interaction, the coupled-cluster singles 
and doubles (CCSD) calculations underbinded substan- 
tially the ground states of ■^~^"He compared to experi- 
ment. However, they correctly predicted ^He and ^He 
unstable with respect to neutron emission. The width of 
the ^He ground state resonance, 0.26 MeV is quite close 
to that calculated here. 

Finally, we note that the NCSMC g.s. resonance en- 
ergy, 0.71 MeV, is lower but still compatible with the 
extrapolated NCSM value of 0.98(29) MeV (see Tables U 
andHni). 



quence, its convergence properties are superior to either 
NCSM or NCSM/RGM. 

We demonstrated the potential of the NCSMC in cal- 
culations of ^He resonances. Starting from a realistic soft 
SRG-N'^LO NN potential that describes accurately two- 
nucleon properties and, with the choice of A = 2.02 fm~^ 
for the SRG evolution parameter, also predicts "^H and 
^He binding energies close to experiment, we calculated 
^He and ^He eigenstates in the NCSM and used them 
as input to the coupled-channel NCSMC equations. We 
found the ^He g.s. energy in very good agreement with 
experiment. The results for the 3/2~ g.s. resonance as 
well as for the well-established 5/2 ~ resonance of ^He are 
in reasonable agreement with experiment. Our results 
for the controversial l/2~ resonance are in fair agree- 
ment with the neutron pick-up and proton-removal re- 
actions experiments of Refs. [H, |49|. Our calculations 
definitely do not support the hypothesis of a low lying 
{Er^\ MeV) narrow (F < 1 MeV) 1/2" resonance. We 
also predict two broad currently unobserved ^^3/2 reso- 
nances at about 3.7 MeV and 6.5 MeV, respectively. 

The NCSMC calculations do not involve any ad- 
justable parameter except for those used in the construc- 
tion of the input iViV (or three-nucleon) potentials. Com- 
putations depend on the size of the HO basis, the HO 
frequency, and the number of eigenstates included in the 
model space. We investigate the convergence behavior 
of the approach with respects to these expansions. Due 
to the over-completeness of the NCSMC basis, the con- 
vergence rate is superior to that achievable with either 
NCSM or NCSM/RGM. The advantages of the NCSMC 
are expected to become even more evident in calculations 
with composite projectiles (such as deuteron, '^H, or '^He) 
that require the use of a large number of pseudostates in 
the NCSM/RGM (or other cluster-based approaches) to 
account for virtual breakup effects. The contribution of 
the pseudostates is expected to be suppressed in the NC- 
SMC approach. Extension of the NCSMC formalism to 
the case of composite projectiles, the inclusion of three- 
nucleon interactions, and the coupling of three-body clus- 
ters are under way. 



IV. CONCLUSIONS AND OUTLOOK 
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Appendix A 

In this appendix we briefly outline the explicit steps for the derivation of the orthogonalized cluster form factors of 
Eq. (|25|) and ([26l) and provide their algebraic expressions. 

The orthogonalized cluster form factor in r-space representation of Eq. (j25p reads 

'M =E / dr'r'\A\rT\A.,^i:j) M~i{r' ,r) 

= E^«^w E ^A^J''T\X'^i7')Ki-^un (Ai) 

= E -^»^('') 9\un , (A2) 

where the orthogonalized cluster form factor in the model-space is given by the model-space non-orthogonalized cluster 
form factor times the model-space norm kernel: 

-gx.n = E {A\.rT\A.'<^iZ)K'l,.n= E 9Xu'n'M;'i,.n- (A3) 

At the same time, the translational- invariant non-orthogonalized cluster form factors in the model space, gx^n, can 
be conveniently derived starting from the Slater-determinant (SD) NCSM eigenstates, 

\AXrT)sD = \AXrT) vooiR^c^l) , (A4) 

and the SD channel states 

\'^i7)sD ^\{\A^aaJ^-T,)sD\aa2irT,))^'^'>YeiRi±)\^' ^^i?„,(41.), (A5) 



and removing the spurious motion of the center of mass. Here, the cm. coordinates of Eqs. (jA4l) and (|A5p are given 
by 

A A 

i?(^1^4fE^- ^cl.-^ E (A6) 

and <y9oo(^c^.) is the HO wave function Rao{Ri'^.)Yoo{Ri'^,). The resulting expression for the non-orthogonalized 
cluster form factor in the singie-nucleon projectile (a = 1) basis is: 

gx^n = {AXrT\A,^ll^) 

^ {nmO, f |00n£, {A\J-T\A.<^fJ} sd 



(n^00,^|00 



SD 



(A7) 



The Moshinsky brackets {n£00,£\00nL£) allows us to transform from the SD to the Jacobi-coordinate states. This 
expression was first derived in Ref. '36l| where further details on the derivation can be found. 
The orthogonalized coupling form factor in r-space representation of Eq. (j26p reads 

hxAr)=Y. J dr'r'\AXrT\H\Ar.>^i:j)M;,i{r',r) 

v' 

= E ^niir) E (.A\rT\il\A.,^iZ) Ar;|, -l- R^_^^r,{r)i,A\rT\AM7^Jkn^aATrei\n^a. + 1^) 

neP i/'n'eP 

= E Rni{r)h\pn + Rnr^a^ + li{r) {nmax£\Trel\nmax + H) gXvn^^^ , (A8) 
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where 

hxun - E {A\rT\H\A.,<i>iZ)N-i^,,,^ h>^'^'n'K'i,.n (A9) 

is the orthogonaUzed couphng form factor in the model space. In deriving the above expression, one has to pay 
attention in taking into account the contribution of transitions to basis states outside of the model space brought 
about by the relative kinetic-energy operator. The model-space non-orthogonalized coupling form factor h\^n can be 
derived in a similar fashion as Eq. (jASp . and is given by: 

hx,n = {A\rT\H\AM7) 

= ^ ^ gXv'n' {n'l'\frel\nl) + gXi^n 
f' n'eP 

+ I V(-i)^i+''-j J ^1 i 

{nmO,£\OOn£,£)_^ ' \ i J j j 2V2 

J'T' (nlj)aba 

■X {{nala3a\,nihjb\)J'T'\V\ inj,j, i , n£j \)J'T') 

X SD{A\rTmal^^^^^,/^^^^^^^^^^^^^ (AlO) 

We note that the point-Coulomb contribution introduced in Eq. (|lll) is omitted in the above expressions for simplicity. 
It is zero in the present application to ^He. Finally, Ei, is t he sum of the eigenergies of the two clusters. 



[1] B. S. Pudliner, V. R. Pandharipande, J. Carlson, S. C. 
Pieper and R. B. Wiringa, Phys. Rev. C 56, 1720 (1997); 
R. B. Wiringa, S. C. Pieper, J. Carlson and V. R. Pand- 
haripande, Phys. Rev. C 62, 014001 (2000); S. C. Pieper 
and R. B. Wiringa, Ann. Rev. Nucl. Part. Sci. 51 53, 
(2001); S. C. Pieper, K. Varga and R. B. Wiringa, Phys. 
Rev. C 66, 044310 (2002). 

[2] P. Navratil, V. G. Gueorguiev, J. P. Vary, W. E. Ormand 
and A. Nogga, Phys. Rev. Lett. 99, 042501 (2007). 

[3] G. Hagen, T. Papenbrock, D. J. Dean and M. Hjorth- 
Jensen, Phys. Rev. Lett. 101, 092502 (2008). 

[4] G. Hagen, M. Hjorth- Jensen, G. R. Jansen, R. Machleidt 
and T. Papenbrock, Phys. Rev. Lett. 108, 242501 (2012). 

[5] R. Roth, S. Binder, K. Vobig, A. Caici, J. Langhammer 
and P. Navratil, Phys. Rev. Lett. 109, 052501 (2012). 

[6] E. Epelbaum, H. Krebs, D. Lee and U. -G. Meissner, 
Phys. Rev. Lett. 106, 192501 (2011). 

[7] K. M. NoUett, S. C. Pieper, R. B. Wiringa, J. Carlson 
and G. M. Hale, Phys. Rev. Lett. 99, 022502 (2007). 

[8] S. Quaglioni and P. Navratil, Phys. Rev. Lett. 101, 
092501 (2008). 

[9] G. Hagen and N. Michel, Phys. Rev. C 86, 021602 (2012). 
[10] K. Wildermuth and Y. C. Tang, A unified theory of the 

nucleus, (Vieweg, Braunschweig, 1977). 
[11] Y. C. Tang, M. LeMere and D. R. Thompson, Phys. Rep. 
47, 167 (1978). 



[12] T. Fliessbach and H. Walliser, Nucl. Phys. A377, 84 
(1982). 

[13] K. Langanke and H. Friedrich, Advances in Nuclear 
Phystcs, edited by J. W. Negele and E. Vogt (Plenum, 
New York, 1986). 

[14] R. G. Lovas, R. J. Liotta, A. Insolia, K. Varga and D. S. 
Delion, Phys. Rep. 294, 265 (1998). 

[15] H. M. Hofmann and G. M. Hale, Phys. Rev. C 77, 044002 

(2008) . 

[16] P. Navratil, J. P. Vary and B. R. Barrett, Phys. Rev. 

Lett. 84, 5728 (2000); Phys. Rev. C 62, 054311 (2000). 
[17] S. Quaglioni and P. Navratil, Phys. Rev. C 79, 044606 

(2009) . 

[18] P. Navratil, R. Roth and S. Quaglioni, Phys. Rev. C 82, 
034609 (2010). 

[19] P. Navratil, R. Roth, and S. Quaglioni, Phys. Lett. B 

704, 379 (2011). 
[20] P. Navratil and S. Quaglioni, Phys. Rev. C 83, 044609 

(2011) . 

[21] P. Navratil and S. Quaglioni, Phys. Rev. Lett. 108, 
042503 (2012). 

[22] P. Descouvemont, E. Tursunov, and D. Baye, Nucl. Phys. 

A765, 370 (2006). 
[23] S. Quaglioni, P. Navratil, G. Hupin, J. Langhammer, 

C. Romero-Redondo, and R. Roth, Few-Body Syst. 53 

(2012) , DOI 10.1007/s00601-012-0505-0. 



16 



S. K. Bogncr, R. J. Furnstahl and R. J. Perry, Phys. Rev. [40] 
C 75, 061001 (2007). 

R. Roth, S. Reinhardt and H. Hergert, Phys. Rev. C 77, [41 

064003 (2008). [42 

R. Roth, T. Neff, and H. Feldmeier, Prog. Part. Nucl. [43 

Phys. 65, 50 (2010). [44 

S.K. Bogner, R. J. Furnstahl and A. Schwenk, Prog. Part. [^45 

Nucl. Phys. 65, 94 (2010). [46 

D. R. Entcm and R. Machlcidt, Phys. Rev. C 68, 041001 [47 
(2003). 

R. Machleidt and D. R. Entem, Phys. Rept. 503, 1 [48 

(2011). [49 

S. Quaglioni, P. Navratil, R. Roth and W. Horiuchi, J. [50 

Phys.: Conf. Ser. 402, 012037 (2012), DOT 10.1088/1742- [5l' 

6596/402/1/012037. [52 

P. Navratil, S. Quaglioni, I. Stetcu and B. R. Barrett, J. [53 

Phys. G 36, 083101 (2009). [54 

S. Baroni, P. Navratil and S. Quaglioni, Phys. Rev. Lett. [55 
110, 022505 (2013). 

P. Navratil and W. E. Ormand, Phys. Rev. C 68, 034305 [56 

(2003). [57 

E. D. Jurgcnson, P. Navratil and R. J. Furnstahl, Phys. [58 
Rev. Lett. 103, 082501 (2009). 

E. D. Jurgenson, P. Navratil and R. J. Furnstahl, Phys. 

Rev. C 83, 034301 (2011). [59 

P. Navratil, Phys. Rev. C 70, 054324 (2004). [60 
M. Hcssc, J.-M. Sparcnbcrg, F. Van Raemdonck, and D. 

Baye, Nucl. Phys. A640, 37 (1998); M. Hesse, J. Roland, [61 
and D. Baye, Nucl. Phys. A709, 184 (2002). 

R. H. Stokes and P. G. Young, Phys. Rev. Lett. 18, 611 [62 
(1967); Phys. Rev. 178, 2024 (1969). 
Z. X. Cao et al, Phys. Lett. B 707, 46 (2012). 



A. A. Korsheninnikov et. al, Phys. Rev. Lett. 82, 3581 
(1999). 

K. Markcnroth et al, Nucl. Phys. A 679, 462 (2001). 
M. Meister et al, Phys. Rev. Lett. 88, 102501 (2002). 

F. Skaza et al, Phys. Rev. C 73, 044301 (2006). 
N. Ryezayeva et al, Phys. Lett. B 639, 623 (2006). 
V. Lapoux et al, J. Phys. Conf. Ser. 49, 161 (2006). 

H. G. Bohlcn et al, Phys. Rev. C 64, 024312 (2001). 

G. V. Rogachev et al, Phys. Rev. Lett. 92, 232502 
(2004). 

A. H. Wuosmaa et al, Phys. Rev. C 72, 061301 (2005). 
A. H. Wuosmaa et al, Phys. Rev. C 78, 041302 (2008). 
D. H. Denby et al, Phys. Rev. C 78, (2008) 044303. 
Yu. Aksyutina et. al. Physics Letters B 679, 191 (2009). 
P. Boutachkov et al, Phys. Rev. Lett. 95, 132502 (2005). 

F. Beck et al, Phys. Lett. B 645, 128 (2007). 

S. Cohen and D. Kurath, Nucl. Phys. A 73, 1 (1965). 

L Brida, Steven C. Pieper, and R. B. Wiringa, Phys. Rev. 

C 84, 024319 (2011). 

R. B. Wiringa, private communication. 

L J. Thompson, private communication. 

I. J. Thompson and F. Nunes, Nuclear Reactions for As- 
trophysics, Cambridge University Press, Cambridge, UK 
(2009) p. 301. 

D. R. Tilley et al, Nucl. Phys. A 708, 3 (2002). 

T. Myo, K. Kato, and K. Ikcda, Phys. Rev. C 76, 054309 

(2007). 

J. Wurzer and H. M. Hofmann, Phys. Rev. C 55, 688 
(1997). 

G. Hagen, D. J. Dean, M. Hjorth-Jensen, and T. Papen- 
brock, Phys. Lett. B 656, 169 (2007). 



